First principle theory of correlated transport through nano-j unctions 



O 

o 



X3 
O 

o 



> 

(N 
(N 
On 
O 

O 



O 

o 



X 
J3 



A. Ferretti/ A. Calzolari,i R. Di Felice/ F. Manghi/ M. J. Caldas,^ M. Buongiorno Nardelli,^- and E. Molinari^ 

INFM National Center on nano Structures and bioSystems at Surfaces (S^) and Dipartimento di Fisica, 
Universita di Modena e Reggio Emilia, ^1100 Modena, Italy 
^Instituto de Fisica, Universidade de Sao Paulo, 05508-900 Sao Paulo, Brazil 
^Department of Physics, North Carolina State University, Raleigh, NC 27695, USA 
^CCS-CSM, Oak Rtdge National Laboratory, Oak Ridge, TN 37831, USA 
(Dated: February 2, 2008) 

We report the inclusion of electron-electron correlation in the calculation of transport properties 
within an ab initio scheme. A key step is the reformulation of Landauer's approach in terms of 
an effective transmittance for the interacting electron system. We apply this framework to analyze 
the effect of short range interactions on Pt atomic wires and discuss the coherent and incoherent 
correction to the mean-field approach. 



One of the most pressing problems in nanotechnology 
is the need for recasting all the know-how about meso- 
scopic transport physics into the fully quantum mechan- 
ical limit appropriate for atomic scale phenomena. In 
the case of electronic and transport properties of atomic 
and molecular conductors, we must address at the micro- 
scopic level both the chemical complexity of the conduc- 
tor and the complexity of the interactions between the 
different components of an extended open system. 

The combination of Green's function methods with a 
Density Functional Theory (DFT) description of elec- 
tronic states has become a standard approach to study 
charge transport at the nanoscale However, 
care must be taken in comparing the computed trans- 
port characteristics to experiments . Indeed, some im- 
portant features — such as electron correlations, dissi- 
pation, decoherence and temperature effects — are at 
the moment partly or fully neglected. These deficien- 
cies become more and more crucial when the dimensions 
of at least one part of the system reach a confinement 
situation where for instance electron-electron (e-e) inter- 
actions may become dominant. Recent observations of 
Kondo effect and Coulomb blockade in molecules con- 
nected to external electrodes 0, 0; ^ind in nanotubes 
with magnetic impurities j3| or size confinement [loj . in- 
dicate that correlations do play an important role in the 
mechanism of charge transport in nano-devices. Whereas 
many efforts |^ 0, 0, Q| have been directed to 
study other aspects of the transport problem, a standard 
approach to include correlation effects does not yet exist. 
One may expect that some of the difficulties in the in- 
terpretation of transport experiments on simple atomic 
chains and individual molecules could be ascribed to the 
neglect or mistreatment of correlations. 

In this Letter we develop a new method for the ab 
initio computation of quantum transport in the strong 
correlation regime, and then apply it to specifically ad- 
dress the effect of e-e interactions on electronic transport 
through atomic-scale conductors. Following closely Meir 
and Wingreen , we derive a Landauer-like expression 
for the current through a correlated conductor between 



uncorrelated leads. Short range e-e interactions in the 
conductor are described thr oug h the Three-Body Scat- 
tering (3BS) formalism [lil flTj. The method is imple- 
mented through use of "maximally-localized" Wannier 
functions H [ll (MLWF). A particularly well-adapted 
system for gauging the relevance of e-e correlations is 
a late-transition-metal break-junction: in such a config- 
uration, e-e effects are negligible in the bulk but, as a 
consequence of dimensionality, may acquire relevance in 
the junction region. We apply our method to simulate 
a model Platinum break-junction of varying length. Our 
findings show a large suppression of the transmittance 
that we ascribe to the inclusion of e-e elastic decoher- 
ence in the simulation. Moreover, the strong conduc- 
tance reduction, also for small wire lengths, suggests that 
correlation cannot be neglected when studying transport 
properties of systems with localized electrons. 

The system is modelled as three different regions, the 
left (L) and right (R) electrodes and a central conductor 
(C) . We express our operators in a localized basis set 0| 
which allows us to write the Hamiltonian and the Green's 
functions of the whole system as 3 x 3 block matrices 
defined on the basis in the L, R, and C regions [H^y, 
Gxy{ijj) where x,y — L,C,R]. 

The Hamiltonian reads: 



J2 Hii'clci' + H,, 



Hii 



h.c. 



(1) 

where c; and cj {di and d|) are the one-electron annihila- 
tion and creation operators in the leads (conductor). In 
the above expression, the first term describes the L and 
R leads, Hint the conductor, the last term the coupling of 
the conductor with the L and R leads. We stress that the 
leads and the coupling Hamiltonian are non-interacting 
and all the e-e interaction is restricted to the conductor. 

From the continuity equation for the steady-state cur- 
rent in the system and using the Keldysh non-equilibrium 
Green's function formalism [i^, the following expression 



for the current [15| is derived: 

y'dc^Tr{[S<-S<]Ac + [ri-rH]G<}. (2) 

Here Ac = i[G^c ~ ^c] spectral function and 

Qr^a,<,> ^j^g (retarded, advanced, lesser, greater) 
Green's functions in the conductor. The trace should be 
taken on the conductor. The interaction between conduc- 
tor and leads is described through the lead self-energies 
(SB's), defined as ^l''''<'> (uj) = Hcx Gl^''^<'> (lu) H^c, 
where x=L,R. Finally, the Tlh terms in Eq. |2Jl are de- 
fined as twice the imaginary part of the retarded lead- 
SE's, i.e. Tl.r = i - J . Note also that 

^L,R = '^fL,R ^L,R and = - Jl.r) Tl,r, where 
Jl and Jr are the Fermi occupation functions for the left 
and right leads. 

While in the non-interacting case the above expression 
brings to the usual Landauer formula , in the presence 
of interaction between electrons this is no longer true 
and further assumptions are needed. We here adopt the 
ansatz proposed by Ng which relates (S^) to 
thus defining the statistics of energy levels, in the 
general out-of-equilibrium interacting case. The starting 
point is the assumption [j^ : 

E<'>H = S<t>HAH, (3) 

where S^iJ^* — S^'^ -I- S^'^ refer to the non-interacting 
case and include only the coupling with the leads, while 
the full S^'^ include also e-e interactions. A(aj), which 
is called A in other formulations is determined by 
the identity E> - E< = S;; - leading to: 

A(.;) ^ [ESc(c^) - Kci^T' i^hiu^) - E^(c.)] . (4) 
Here the interacting SE's take the form: 

E^"(C.) = + ^Ti^) + Ko^rH (5) 

where TT'^orr account for the e-e interactions (while Eq^ 
just drop this last term). Relations I|3I5(I are the key 
to relate Eq. Q to a Landauer-like formula. In fact, 
following Eq. (6) of Ref. and using = GJ;E<G^, 
it is possible to derive G^ = iG^f/iFL + /hFtj] AG^. 
and therefore G^ - G^ = -i G^IFl F^] AG;i;. These 
steps lead to the final expression for the current: 

^[/L-/i?]Tr{FiG^F«AG^}. (6) 

We remark that the e-e correlation plays a twofold role: 
it renormalizes the Green's functions, which should now 
be calculated for the interacting system, and also mod- 
ifies the expression for the current, introducing the cor- 
rective factor A(a;). The quantity traced in Eq. ©, even 
if not a true transmittance across the scattering region 
due to the breakdown of Landauer's theory, still plays 



the same role as regards transport. For this reason we 
refer to it as an effective transmittance and compare it 
to the transmittance of the non-interacting case. In par- 
ticular, since the imaginary part of the e-e self-energy 
exactly vanishes at the Fermi energy {Ep), the correc- 
tion A{uj — Ep) becomes the identity operator and the 
effective transmittance computed at Ep is proportional 
to the conductance at zero temperature, provided that 
the SE's are calculated for the interacting system. 

An accurate evaluation of the effective transmittance 
needs to be based on: (i) a good description of the non- 
interacting system, and (ii) the calculation of the e-e 
correlated SE to include many-body effects arising from 
the interaction in the conductor. The first problem is 
solved by exploiting a recent methodology for the ah ini- 
tio calculation of the transmittance in the coherent trans- 
port regime The ground state of the mean- field sys- 
tem is described within the Local Density Approximation 
to DFT, using norm-conserving pseudopotentials and a 
plane- wave basis set jl^. The geometry of the L-C-R 
nano-junction are as in Ref. [l|. To obtain a consistent 
description of the system in a real-space localized basis 
set, the MLWF's are computed both for the leads and 
for the conductor [i^l ■ The details of this transformation 
and its application to electronic transport arc described 
elsewhere |^ 0| . We remark that the basis change from 
Bloch eigenvectors to MLWF's preserves orthonormal- 
ity and completeness in the original Hilbert subspace, 
thereby avoiding typical problems arising very often for 
other localized basis functions. The same features al- 
low us to use a minimal basis set for computing trans- 
port properties, while employing the system independent 
plane-wave basis set for the DFT calculation. 

To discuss the inclusion of e-e correlations we first need 
to define the interaction hamiltonian in Eq. . We focus 
on the short range e-e interaction for two main reasons. 
On one hand, this regime is characterized by strong de- 
viations from the non-interacting behavior, for instance 
in terms of quasiparticlc lifetimes (which include elastic 
decoherence in the transport formalism). On the other 
hand, it allows us to adopt an Anderson-like form 
of the interaction which is suitable for the localization in 
the conductor region only, as required by our approach. 

In this work the e-e self-energy E^^^^ is computed using 
a non-perturbative approach based on an effective Ander- 
son hamiltonian, whose U Coulomb integrals could be ei- 
ther calculated ah initio 26] or used as adjustable input 
parameters. It is solved by means of a configuration in- 
teraction scheme with up to three bodies (3BS) added to 
the non-interacting Fermi sea: two (one) electron(s) and 
one (two) hole(s) for conduction (valence) states. This 
method has been successfully applied to describe photoe- 
mission experiments on strongly correlated systems |lZi ■ 
The BBS self-energy is formally given 'l7| as a sum over 
projectors onto atomic states (those with non-negligible 
U integrals) and thus can be properly localized in the 



conductor. 

We now come to the model system: a Pt atomic wire of 
varying length, where the correlation is switched-on only 
on a finite number Nc of atoms. It is worth noting that 
late bd transition-metal atomic wires have been demon- 
strated in break-junction experiments Moreover, in 
the case of Pt, correlation effects are expected to strongly 
increase their importance in passing from the bulk leads 
to the confined wire as an effect of the reduced dimen- 
sionality. Therefore, our model system, while neglecting 
the full complexity of the wire-lead interfaces, retains the 
basic geometric construct of a correlated wire between 
uncorrelated leads in a Pt break-junction. Hence, it can 
be considered suitable to emphasize the effects of cor- 
relation and decoherence in nano-junctions, whereas the 
effects of dimensionality at the junctions (3D instead of 
ID leads) is already known 27]. 

The setup is described in Fig. ^a). For region C we 
consider a supercell containing 11 atoms with an inter- 
atomic separation of 3.0 A, reasonable for experiments 
under stretching (such as break-junctions ^27]). The 
same geometry is adopted for the leads, which are mod- 
elled as semi-infinite wires treated at the single-particle 
level. Since the sd-bands of the Pt chain [Fig. Q^b)] 
form an isolated subset jl8l |. we can extract a manifold 
of MLWF's which span the same subspace. As an ex- 
ample, a particular wave function well localized within 
two bond- lengths is shown in the inset of Fig. ^c) . This 
MLWF subspace allows us to reproduce with good ac- 
curacy states more than 2 eV above the Fermi energy, 
enough to describe transport properties in this system 
where p-orbitals are expected to play a negligible role. 
Figure ^c) shows the computed transmittance for the 
non-interacting Pt chain, which counts, for each en- 
ergy, the number of channels (bands) available for charge 
transport, leading to a zero temperature conductance of 
3Go (Go = 2eV/i). 

In Fig.|2Ia,b) we report the spectral function projected 
on the interacting region and the effective transmittance 
for a chain with three correlated atoms (U = 2.0 eV). 
The many-body spectral function shows a splitting of 
the d-bands and a slight upward energy shift, which is 
consistent with our short range interaction picture based 
on the physics of the Anderson hamiltonian. The trans- 
mittance is strongly suppressed by the inclusion of cor- 
relation [FigJ2b)], which is particularly effective in the 
hole region [28j . This can be understood considering that 
on-site correlation arises from the strongly localized Pt 
d-orbitals, that are largely occupied and produce major 
features below Ep- 

The introduction of e-e interactions leads to quasi par- 
ticles characterized by energy and lifetime (finite broad- 
ening of their spectral features) formally accounted for 
by the hermitian {H) and the antihermitian {A) parts 
of the SE operator. Note that in more common correla- 
tion treatments, such as LDA-I-U, the SE in hermitean 
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FIG. 1: (color online), (a) System layout divided in the leads 
and conductor regions. Open (full) circles are uncorrelated 
(correlated) Pt atoms, (b) Band structure of the Pt infinite 
chain. A is the edge of the one-dimensional Brillouin zone, (c) 
Transmittance of the non-interacting conductor region. The 
inset shows the isosurface plot of a computed MLWF. The 
Fermi energy is set to zero. 

and lifetimes are consequently neglected. Although these 
components must obey analytic constraints |r^ . we an- 
alyze them separately to highlight their very effect on 
transport. First of all we divide the SE operator in H 
and A contributions, Yicorr = + "^A where: 

^HA^) = \[n^)±^\^)\ ■ (7) 

In Fig. |2tc,d) and (e,f) we show our results when us- 
ing Y,H or Syi, respectively, instead of the full correlation 
SE from 3BS. By definition Y,a vanishes at Ep, thus 
the only contribution to the zero temperature conduc- 
tance comes from the H part of E. Moreover, at en- 
ergies different from Ep the major quenching factor on 
the transmittance is due to the A part of S, while its 
effect on the spectral function is just a slight broaden- 
ing of the main peaks. Note that the decrease of the 
effective transmittance due to Y^h is related to energy 
misalignment of channels between the conductor and the 
leads, while that due to Y,a accounts for e-e scattering. 
Our results thus indicate that e-e scattering plays a fun- 
damental role in suppressing effective transmittance for 
systems with strong short range correlation. This effect 
should be hidden in the low temperature conductance but 
evident in the I — V characteristics [Eq. a suggestion 
for further experimental analysis. 

Fig.Elreports the computed conductance and effective 
transmittance for variable number of correlated atoms 
Nc- Figure|3Ia) shows the effective transmittance curves 
for various Nq values and demonstrates that the effec- 
tive transmittance decreases with increasing Nc- This 
behavior is consistent with the presence of a scattering 
mechanism, according to which one expects a vanishing 
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FIG. 2: (color online). Curves: Spectral function projected 
on the interacting region (left) and effective transmittance 
(right) for the case of three correlated atoms in the Pt chain. 
Shaded areas: mean-field reference results. Ef is set to zero, 
(a) and (b) show the computational results obtained with 
the total correlation SE Ecor-r = ^3bs- The (c,d) [e,f] plots 
are obtained using only the hermitian [antihermitian] SE E^f 
Pa]. 
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FIG. 3: (color online), (a) Effective transmittance against 
energy for variable number of correlated atoms in the chain: 
A is the reference mean-field bulk transmittance. A^c=l,3,5,7 
in B,C,D,E. (b) Conductance as function of the number of cor- 
related atoms in the Pt chain. Saturation is observed between 
Nc=6 and 7. 

effective transmittance at w ^ Ep in the limit of in- 
finitely long correlated wires. Figure Olb) reveals also 
a rapid decrease of the conductance at small values of 
Nc- The last result shows firstly that the conductance is 
renormalized as well as the transmittance, although the 
imaginary part of the SE vanishes at Ep. Secondly, it 
suggests that the effect of on-site correlation should be 
experimentally measurable also for short atomic chains 
(such as those produced by break-junctions). 

In conclusion, we derived a generalized Landauer-like 
expression for the current [Eq. jnj] and the conductance 
in the presence of e-e interactions. The formalism is suit- 
able for a fully ah initio implementation that we realized 
using a basis set of maximally-localized Wannier func- 



tions for the Green's functions and the 3BS formalism 
for the e-e self-energy. We applied the method to a finite 
Pt wire and found a renormalization of the conductance 
and a strong quenching of the transmittance as a conse- 
quence of the e-e scattering. Our results suggest that the 
inclusion of electron correlation for systems with strong 
short-range interactions is essential for an accurate de- 
scription of current and conductance. 
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